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The vacuum average value of the stress-energy tensor of a massive scalar field with non-minimal 
coupling ^ to the curvature on the short-throat flat-space wormhole background is calculated. The 
final analysis is made numerically. It was shown that the energy-momentum tensor does not vi- 
olate the null energy condition near the throat. Therefore, the vacuum polarization cannot self- 
consistently support the wormhole. 
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I. INTRODUCTION 



Traversable wormholes have been introduced into physical realm by Morris and Thorne [l| in 1988. At the 
same time they realized that a matter threading a wormhole's throat should possess exotic properties, namely 
^ it should have a negative pressure and violate the null energy condition (NEC). Later on, this result was 
^ generalized for any traversable wormholes, both static and non-static 2]. As is known, the classical matter does 
5-H , satisfy the usual energy conditions, hence traversable wormholes cannot arise as solutions of classical relativity 
I and matter. Already in 1989 Morris, Thorne, and Yurtsever Q supposed that quantized fields could play a role 

[ of the exotic matter maintaining wormholes. Their reasons were founded on the important fact stating that 
' quantum field theory may have states with negative energy density, thus violating the weak energy condition 
> ; 13 (see, also, 0). 

\! . In the absence of the complete theory of quantum gravity, the semi-classical approach gives the more natural 
way to include quantized fields in the theory of gravity. Various wormhole solutions in semi-classical gravity 
have been considered in the literature. For instance, semi-classical wormholes were found in the framework 
of the Frolov-Zelnikov approximation for (T^i/)'^''" [6]. Analytical approximations of the stress-energy tensor 
, of quantized fields in static and spherically symmetric wormhole space-times were also explored in Refs. 
[^. Some arguments in favor of the possibility of existence of semi-classical wormholes have been given by 
Khatsymovsky llOll. However, the first self-consistent wormhole solution coupled to a quantum scalar field was 
obtained in Ref. |11| . The ground state of a massive scalar field with a non-conformal coupling on a short-throat 
flat-space wormhole background was computed in Ref. p^ . by using a zeta renormalization approach. The 
latter wormhole model, which was further used in the context of the Casimir effect [l^, was constructed by 
^ ■ excising spherical regions from two identical copies of Minkowski space-time, and finally surgically grafting the 
5^ I boundaries (A more realistic geometry was considered in Ref. [l^)- In a series of works lB"l3 various aspects 



of the graviton one loop contribution to a classical energy in a wormhole background have been analyzed. The 
latter contribution was evaluated through a variational approach with Gaussian trial wave functionals, and 
the divergences were treated with a zeta function regularization. In particular, the finite one loop energy was 
considered as a self-consistent source for a traversable wormhole. 

Note that up to now no one has succeeded in exact calculations of vacuum expectation values of the stress- 
energy tensor of quantized fields on the wormhole background. The reason for this state of affairs consists in 
considerable mathematical difficulties which one faces with trying to quantize a physical field on the wormhole 
background. To overcome these difficulties, in this work we will consider a simple model of wormhole space-time 
given in [l^: the short-throat flat-space wormhole. The model represents two identical copies of Minkowski 
space with excised from each copy spherical regions, and with boundaries of those regions are to be identified. 
The space-time of this model is everywhere flat except a two-dimensional singular spherical surface. Due to 
this fact it turns out to be possible to construct the complete set of wave modes of the massive scalar field and 
calculate the stress energy tensor. 
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The organization of the paper is as follows. In Sec. |TT]we describe a space-time of wormhole in the short-throat 
flat-space approximation. In Sec. IIIII we analyze the solution of equation of motion for massive scalar field and 
obtain Euclidean Green function. Sec. IIVI devoted to analysis of the vacuum expectation value of the square of 
field from analytical and numerical point of views. In Sec. |V] we obtain close formulas for numerical analysis 
the stress-energy tensor and in Sec. IVII discuss our results. 

We use units h ^ c = G = 1. The signature of the space-time, the sign of the Riemann and Ricci tensors, is 
the same as in the book by Hawking and Ellis [l9| . 



II. A TRAVERSABLE WORMHOLE: THE SHORT-THROAT FLAT-SPACE APPROXIMATION 

In this section we briefly consider a simple model of a traversable wormhole (see Ref. 0]). Assume that the 
throat of the wormhole is very short, and that curvature in the regions outside the mouth of the wormhole is 
relatively weak. An idealized model of such a wormhole can be constructed in the following manner: Consider 
two copies of Minkowski space, A4+ and Ai-, with the spherical coordinates {t,r±,9±,(p±) [Notice: A4+ and 
M- have a common time coordinate t. One may interpret this fact as the identification t+ o t^.]; excise 
from each copy the spherical region r± < a, where a is a radius of sphere; and then identify the boundaries of 
those regions: (<, a, <y9_|_) O (t, a, The Riemann tensor for this model is identically zero everywhere 

except at the wormhole mouths where the identification procedure takes place. Generically, there will be an 
infinitesimally thin layer of exotic matter present at the mouth of the wormhole. 

Such an idealized geometry can be described by the following metric 

ds^ = -dt^ +dp^ +r^{p) {d0^ +3111^6 difi^), (1) 

where \p\ is a proper radial distance, — oo < p < oo, and the shape function r(p) is 

r{p) = \p\+a. (2) 

It is easily to see that in two regions 72.+ : p > and TZ- : p < separately, one can introduce a new radial 
coordinate r± = ±p + a and rewrite the metric ([1]) in the usual spherical coordinates: 

ds^ = -dt^ + drl + r\ (de^ + sin^ Odif^), (3) 

This form of the metric explicitly indicates that the regions TZ+ : p > and TZ- : p < are fiat. However, note 
that such the change of coordinates r = |p | + a is not global, because it is ill defined at the throat p = 0. Hence, 
as was expected, the space-time is curved at the wormhole throat. To illustrate this we explicitly calculate the 
scalar curvature R{p) in the metric ([1]): 

R{p) - -8a-'6{p). (4) 

HI. GREEN FUNCTION 

Let us discuss a massive scalar field (f> with non-minimal coupling to curvature. The scalar field equation of 
motion has the following form 

{a-^R^m'^)(j)^0. (5) 
A corresponding equation for the Euclidean Green function reads 

(a-m^-ei?)G£(x;i) = -^^p^, (6) 

V9 

where t = —it is the Euclidean time. Due to spherical symmetry one can represent Ge{x; x) as follows 

oo -i-i „ , 

G{x;S:) = J] YC^imM^) / ^e^^^^ g^,{p, p). (7) 

1=0 m=-l •' 

Here gi^i (p, p) is a radial Green function obeying the equation 

9ui + —9ui - ( A + + 1 g^i = — . (8) 
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where = w^ + m^ and a prime means the derivative with respect to p. A solution of Eq. ([8]) can be represented 
as 



9ojiip,p) ^ Q{p - p)<i>+{p)<i>-{p) + Q{p ~ p)'t>+{p)<i>-ip), 

where (j)± (p) are two independent solutions of the radial equation of motion: 



2r' 
b" + —<tJ 
r 



A2 + fcll+ei?l0 = O. 



(9) 



(10) 



In general, Eq. (jlOp cannot be solved in an explicit form. However, this is possible in the particular case 
''(p) — \p\ + this case R = —8S{p)/a and Eq. (|T0|) yields 



_ ^Sip) U = 



(11) 



Note that Eq. ([TT|) is an ordinary second-order differential equation with a delta-like coefficient, S{p). A common 
treatment of such equations is to solve them separately in regions with p 7^ and then match the obtained 
solutions and their first derivatives at p = 0. In the regions p 7^ a general solution of Eq. ([TT|) reads 

(12) 



V zr 

where Ci and C2 are constants of integration and v ~ I + \/2. Integrating Eq. (fTT|) around p = gives the 
following matching conditions at the throat 



4>{+Q) - 4>{-0) = 0, 

(/.'( + O)-0'(-O) : 



</>(+o). 



(13) 



Using the formula ((T^) and the relations ([T^ one may define constants of integration Ci and C2 and then 
construct the radial Green function g{p,p) given by Eq. ([5]). As the result, one obtains (see for details Ref. 



gLoi{p,p) = gt\{p-,p) 

if p are p have the same signs, and 

gcuiip,p) = 
if p and p have different signs. Here 



Xa{I,Kl + r^K,) + (8f - l)hK, 



2XaK^Kl + (Se - l)K^ 



K,{Xr)K,{Xf) 



Xa 



2XaK,Kl + (8e - l)K^ 

K,{Xr)I,{Xf) 



9t\{P,P) = 



(14) 



(15) 



(16) 



is the radial Green function of the Minkowski spacetime. 



IV. VACUUM POLARIZATION 



A vacuum of any quantized physical field is polarized in curved space-time. A vacuum polarization (0^) of the 
massive scalar field in a wormhole space-time has been discussed in Refs. 0, Q in the WKB approximation. 
In this section we represent exact calculations for {(jp') in a wormhole space-time with the metric ([T]). The 
renormalized expression for the vacuum polarization is defined as follows 



lim [Ge{x] x) - Gds{x\ x)] , 



(17) 



where Gds{x',x) are the well-known DeWitt-Schwinger counterterms |2Q |: 



Gds{x]x) 



Ai/2 



87r2 



ao 
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rn'^L ( 1 + —rn^a + 



-TTi^fT 
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a2U 
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+ 



(18) 
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where a is half the squared distance between the points x and x along the shortest geodesic connecting them, 
L = 7 + ^ In !2-J£l and 7 is Euler's constant. In the above expression we have to take all terms which survive 
after taking two derivatives and going to the coincidence limit. Here are the heat kernel coefficients (see, for 
example [2l|). 

The spacetime under consideration (JlJ possesses the singular curvature (j4]) and therefore we have a spacetime 
with the singular surface, p = 0, with codimension one. In this case we have to use expressions for heat kernel 
coefficients obtained in Ref. (2^ (see also Ref. [llj). The main result of the paper [22| is that the heat 
kernel coefficients in this case may be represented as a sum of standard expressions for heat kernel coefficients 
calculated without a singular surface and surface terms: 



a„(M,E) =a„(Af\E)+a^(S), 



(19) 



where M is a manifold, S is a singular surface and Oq (S) — 0. The spacetime given by the metric ([T]) is flat 
and therefore ao(M \ S) = 1, a„>i(Af \ S) = 0. All surface terms a^(I]) are localized at the singular surface. 
For example 
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a|(S) = |?^^<5(p). 
Therefore for p 7^ all singular contributions are zero and we find 

1 _2 



Gds{x]x) 



87r2cr 



ra 
8^ 



1 , m^lal 



m, 

16^ 



(20a) 
(20b) 

(21) 



In fact, this expression for Gds{x]x) does not contain any curvature terms and therefore coincides with that 
of Minkowski space-time. Because of this fact the renormalization procedure reduces to discarding the term 
5*1 (P:P) in Eq. Now, by using the formulas Q, ^ and ([21]), we obtain 



1 C°° °° 



Xa{I,Kl + IlK,) + (8e - l)I,K, 



2XaK^Kl + (Se - l)K^ 



Kl{\r). 



(22) 



For p = —(1 + 0/2 = 4f — 1 we observe singularity in integrand for Z = for = {4^ — 1)^ — m?a?'. It is 
possible for ^ > 1/4. For ^ < 1/4 there is no singularity. We will assume ^ < 1/4. 

Another singularity is for p = 0. Indeed for p = the uniform expansion of integrand gives 



(0'(P)> 



1 p ^ n 

•^0 1=0 ^ 



'(1 + C - 1") - —t\i + c - t^f + o{v-^) 



(23) 



with t — l/Vl + z^. The first two terms of the series are divergent. 

The result of numerical analysis of {(fP') given by Eq. is shown in Figs. [T]and[5] 




FIG. 1: Plots of {(f?) a? for ^ = 0, |, |, 0.2, 0.4 from bottom to up. The value of the field mass m is fixed: m — 1. The 
(</)^) behaviour is shown on the figure b for small distances at the throat. All lines fall down to minus infinity near the 
throat. 
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FIG. 2: Plots of (0^) for m = 0, 1, 2 from bottom to up. The value of $ is fixed: C = 0. 



V. VACUUM STRESS-ENERGY TENSOR (T^^) 
To calculate the stress-energy tensor one may use the standard formula 



where G is the renormalized Euclidean Green function. For above function it is easy to show that 



[GJ 






= 0, 9^ p 




= [G,^'i/'], 




= [G;^V], 







For this reason we may rewrite the expression for EMT in the following form 



(1 - 20G;^.' - 2CG;^. + ( 2^ - 2 ) + K + ^G) + ?i?^.G 



For calculation the EMT we have to take into account that = and P/(l) = l{[ + l)/2. For brevity let 
define operator £ and the function / 



1 



27r2 ./n 2XaK,Kl + (8^ - l)^!:^ 



{•} 



Aa 



The components of the EMT have the following form 



{TP = c\{i- 2or - ^Uf + 2^ - 



{TD = t 



1(1 + ^)^2 2e 



2r2 



2 



where an overdot notes the derivative with respect to p. The function / obeys to the equation 



.. 2 . + 



A' / = 0. 
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VI. NUMERICAL ANALYSIS 

The expressions p8 | - ([3T|) for components of (T^^) are not much suitable for an analytical consideration, 
therefore, we have applied numerical methods for their analysis. In this section we will discuss results of 
numerical analysis. 

For a static spherically symmetric configuration one has (Tj*) = — e, {Tp — p, and (Tg) — {T^) — pt, where 
£ is the energy density, p is the radial pressure, and pt is the transverse pressure. Values of £, p, and pt are 
connected by the conservation law {Tlf)-f^ = 0, which for the Minkowsky metric ([3|) takes the simple form: 

Pt=p+^rp'. (32) 

From here one may easily find pt provided p is found. 




FIG. 3: Plots of pa* for m = 1 and ^ = 0, ^, ^, 0.2, 0.3. The thicker is the line, the grater is the value of ^. 

We compute numerically the energy density e — — (T/) and the radial pressure p — {Tp using Eqs. ([28)) . 
([29| . Results of numerical computations are given in Figs. [SHSl Let us discuss them in details. In Fig. [3] plots 
of p{p) are shown for a fixed value of m and various values of the curvature coupling parameter ^. It is seen that 
p is everywhere negative, and p — > —oo in the limit p — > 0. Thus, a vacuum polarization leads to an infinitely 
negative radial pressure p at the wormhole's throat. 



EC* 




FIG. 4: Plots of ea* for m — 1 and ^ = 0, |, ^,0.2,0.3. The thicker is the line, the grater is the value of ^. The ea* 
behaviour is shown on the figure b for small distances at the throat. All lines fall down to infinity near the throat. 

Plots of energy density e{p) are shown in Fig. |4l It is seen that a qualitative behavior of e depends on ^. 
Provided ^ < 1/8 or ^ > 0.2, the function e{p) reaches a negative minimum at some p, and tends to zero far 
from the throat. In case 1/8 < ^ < 0.2 e is everywhere positive monotonically decreasing function. It is worth 
noting that in both cases the vacuum energy density goes to infinity at the wormhole's throat, i.e., e — )> cx) in 
the limit p (see Fig. |41d). 

It is particularly important for a wormhole geometry to check whether the vacuum stress energy tensor obeys 
the usual energy conditions. In particular, the null energy condition (NEC) reads {Tf^^)k'^k'^ > 0, where fc'^ is 
an arbitrary null vector. In a static spherically symmetric case the NEC reduces to e + p > 0. Graphs for the 
combination e + p are given in Fig. [5] It is seen that the combination e + p behaves similarly to the energy 
density e. In particular, e + p ^ oo in the limit p — > (see Fig. [5}d). Thus, the vacuum stress energy tensor 
(T^i,) does not violate the NEC in the vicinity of the wormhole's throat. 
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FIG. 5: Plots of (e + p)a'^ for m = 1 and ^ = 0, i, i, 0.2, 
(e + p)a'* behaviour is shown on the figure b for small dii 
throat. 
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VII. CONCLUSION 



We have calculated the vacuum polarization (0^) and components of the vacuum stress energy tensor (T^^) 
of the massive scalar field in a wormhole spacetime using the short-throat flat-space approximation for the 
wormhole geometry. The most important result obtained consists in the fact that (Tf^^) does not violate the 
NEC in the vicinity of the wormhole's throat. As a consequence, this implies that the vacuum polarization 
cannot self-consistently support the wormhole. Of course, it is necessary to emphasize that this conclusion has 
been obtained for the very simple model of wormhole. To make more general conclusions one should consider 
more realistic models of wormholes. 
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